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COSMIC GAMMA-RAYS FROM 


METAGALACTIC PROTON-ANTIPROTON ANNIHILATION 

F. W. Stecker* 

ABSTRACT 

The only practical way of determining the existence of cosmological anti- 
matter is by the detection of gamma-radiation produced by the mutual annihila- 
tion of matter and antimatter in the universe. It is shown here that annihilation 
gamma -radiation has a characteristic spectrum which distinguishes it from 
other types of cosmic gamma -radiation. Such spectra are presented here, based 
on the results of recent accelerator experiments on antiproton annihilation and 
also including cosmological distortions of local annihilation-gamma-ray spectra. 
Using the "baryon inhomogeneity" cosmology of Harrison, it may be possible to 
account for both the observed cosmic isotropic X-rays and gamma -rays as due 
to redshifted annihilation-gamma-radiation where the mean value of the product 
of interacting matter and antimatter densities is of the order of UP 18 cm' 6 . 
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COSMIC GAMMA-RAYS FROM 


METAGALACTIC PROTON-ANTIPROTON ANNIHILATION 

F. W. Stecker* 


1. INTRODUCTION 

The existence or nonexistence of antimatter in the universe is a question of 
importance in the fields of cosmology and particle physics (Alfven 1965). 1 Re- 
cently, Harrison has proposed the existence of large amounts of antimatter in 

the universe in order to account for the condensation of matter into galaxies 
2 

(Harrison 1967). The only practical way of determining the existence of such 
antimatter is by the detection of the gamma -radiation produced when antimatter 
and matter annihilate into mesons having gamma-ray -producing decay modes. 

The gamma -ray spectra from such annihilations are presented here. They are 
based on the results of recent accelerator experiments on antiproton annihilation 
and also include cosmological distortions of local annihilation-gamma-ray 
spectra. 

2. THE LOCAL ANNIHILATION- GAMMA -RAY SPECTRUM 

The annihilation gamma -ray spectrum (AGS) for proton-antiproton interactions 
involving kinetic energies less than 286 Me V (the threshold for non-annihilation 


*NRC-NASA Resident Research Associate. 


1 



pion production) is given by the expression 


MM = B J d '’'' f ('')2]f dE . cr =( E . ;v )2]MV f ^(v E J (i) 

s d 

where v is the relative velocity of the proton and antiproton, f(v) is the nor- 
malized distribution function over this velocity, s is an index representing the 
particular type of particle produced in the annihilation, a s ( E s ; v) is the cross 
section times multiplicity of particles of type y and energy E s produced in a 
collision of velocity v, the index d specifies a specific decay mode with a branch- 
ing ratio, R yd , which produces £ yd gamma-rays with a normalized energy dis- 
tribution, f(E y ; E s ). The quantity, B, is defined as the product of interacting 
proton and antiproton densities integrated over a line-of-sight path-length^, i.e.. 


B 



^ n P (^) n p (^) 


( 2 ) 


It can be shown that for energies less than 286 MeV, the production cross 
section, cr g ( E s ; v) can be written as 

ct s( e s ; v ) = °- s £ s V_1 cf s ( E s) (3) 

where f s (E s ) is again a normalized distribution function and a, and Z_ s are ap- 

3 

proximately constant production cross sections and multiplicities (Stecker 1967). 


2 



With the help of Equation (3), Equation (1) reduces to 


mm = ■ (4) 

s ^ d 

Fortunately, it is not necessary to use the general form of Equation (4) in 

order to evaluate the local AGS. Recent accelerator data on proton-antiproton 

annihilations indicate that only about 20% of the gamma -rays produced arise 

through nonpionic meson production (see references and discussion in Stecker 
3 

1967). Furthermore, the largest nonpionic contribution to the AGS is due to the 
/3-meson decay schemes 


P 


± 


77 * + 77 ° 

y + y 


( 5 ) 


(The decay scheme, p° ->77° + 77 0 , is forbidden by the conservation of isospin in 
strong decays.) Gamma-rays from /3-meson decay possess an average energy of 
210 MeV, not much different from the 190 MeV average energy given to gamma- 
rays from the directly produced pions that are an order of magnitude more 
frequent. Other mesons being even less frequently produced than the p -me sons, 
we can conclude that mesons other than pions have a negligible effect on the 
total gamma-ray spectrum from proton-antiproton annihilation. 

We need therefore include in Equation (4) , only the contribution from neutral 
pion decay, thus eliminating the summation over s . The form of the distribution 

, \ 4 

function, f ds (E y ; E s ) is well known (see, for example, Stecker 1966), and there 


3 



is only one important decay mode to consider, viz., n° - 2 y, thus eliminating the 
summation over d. Equation (4) then reduces to 



2B a l 

TT ° 7 7 




( 6 ) 


3. PION PRODUCTION IN PROTON-ANTIPROTON ANNIHILATION 


It now remains only to evaluate the quantities cr. and ^ and the production 
function (E^) in Equation (6). Of the annihilation processes which yield 
neutral pi-mesons, the process 


P + P - 77° (7) 

is, of course, forbidden by conservation of momentum. The process 

p + p -* 77° + 77° (8) 

is also forbidden when it is noted that proton-antiproton annihilations near rest 

occur predominantly from the S states of the proton-antiproton system. It has 

5 

been shown by Lee and Yang (1956) that reaction (8) is then forbidden by the 
conservation of G-conjugation parity (see also Stecker 1967 for a more detailed 

3 

discussion) . Therefore, the extremum gamma -ray energies that one would 
ejqject from proton-antiproton annihilation would result from annihilations of 
the type 


p + p -> 77 + +77 + 77 ® 


(9) 


4 



which produce neutral pi-mesons with a maximum c.m.s. energy of approximately 
870 MeV. Such mesons can produce gamma -rays only within the energy region 

5 MeV < E y < 865 MeV . (10 ) 

(The electromagnetic process, p + p - 2 y may also occur, but is rare enough to be 

4 

neglected (Stecker 1967). 

The product, 2cr n £ , has been determined experimentally to be 

2°V4 - 1-44 x 10' 2 5 cm 2 (11) 

4 

(see references in Stecker 1967). 

6 

The function f n (e^) , was taken from the calculations of Maksimenko (1958), 
based on the statistical theory of multiple particle production. It has been shown 

7 

(Matsuda 1966) that the statistical theory of multiple particle production pre- 
sents excellent agreement with the data in describing annihilations into three or 
more mesons. 

The resulting gamma-ray spectrum, up to 750 MeV is shown in Figure 1. 

8 

Frye and Smith (1966) have calculated the form of this spectrum up to 500 MeV, 
based on recent measurements by the Columbia University group on charged 
pion production in pp annihilation. The excellent agreement between the results 
of Figure 1 and the calculations of Frye and Smith not only serves as a mutual 
check on the calculations, but also supports the conclusion that mesons other 
than pions have a negligible effect on the total gamma-ray spectrum. 


5 



4. THE METAGALACTIC ANNIHILATION-GAMMA-RAY SPECTRUM 


The metagalactic annihilation-gamma -ray spectrum will differ from the 

local spectrum calculated in the previous section due to contributions at non- 

9 

zero redshifts, an effect which has been discussed previously (Stecker 1968). 
Therefore, Equation (6) must be modified to include the cosmological effects of 
curvature, volume and density changes in an expanding universe, time dilation, 
and energy redshift. In a Friedmann -type expanding universe with a Robertson- 
Walker metric the proper relativistic formula for calculating the spectrum 
becomes 


I. (E ) = 2a £ c 

A \ y / 7 T ^7T 


f r [a + *)£,] de 


WZ m a x fy [(1 + Z 

dz n (z) n_ (z) 

J 0 P P (1 + z) 


y- 

3 dz 


where 


(M ■ f 


870 MeV 


dE 


( E .) 




2 

'77 


and z is the cosmological redshifto The curvature factor, di/dz is given by 


di 

dz 


cH " 1 


(1 + z) 2 (1 + 2q 0 z) 


1/2 


(14) 


(Sandage 1961) 10 where c is the speed of light, H 0 is the present value of the 
Hubble parameter and q 0 is the deceleration parameter. We consider the values 
q 0 — 0 and q 0 = 1/2, corresponding to the low density and Einstein-de Sitter 


6 



cosmological models respectively. We also assume that the densities of inter- 


acting matter and antimatter are inversely proportional to the volume of the 

universe and therefore are proportional to (1 + z) 3 » Such would be the case if 

matter and antimatter were effectively separated by a primordal baryon inhomo - 

2 

geneity (Harrison 1967) so that only a small fraction annihilates through peri- 
pheral interactions involving these inhomogeneities . Under this assumption, 
Equation (12) reduces to 


i»(M - 


2B 0 int 




f 


dz(l + z) a f [(l + z)E ] 


(15) 


where 


B 


int 

0 


n int n- int 
p,o 1 p.o 


cH “ 1 


(16) 


nj”* and n_ in g are the present densities of interacting matter and antimatter 
respectively, and the index a is equal to 1 for the low density model and 1/2 for 
the Einstein-de Sitter model of the universe. 

It is of interest to note that since f y (E ) is nonvanishing only within a 
restricted energy range (given by Equation (10)), there exists an energy region 


kJ mc v ^ 

(l + z 1 ^ E v < 5 MeV (17) 

where the integral in Equation (15) is independent of the form of f y ( E ry ) and the 
resultant spectrum has power-law form. In actuality, the power-law form of the 
spectrum is a valid approximation even over a wider energy range where a change 


7 



of the limits of integration over f y (E y ) does not significantly affect the value of 
the integral. In this case, we find 

I. ( E ) = 2B n int a l cKE'< a+1 > ~ E"( a+1 ) 

A \ y ) 0 rr b 7r y y 

where 

1 *865 MeV 

dE y E a y f y (E y ) (1 8) 

5 MeV 

Equation (15) was used to evaluate the form of the metagalactic annihilation- 

gamma-ray spectrum for various values of z max . These spectra are shown in 

Figures 2 and 3 for the Einstein-de Sitter universe and Figures 4 and 5 for a 

low density universe. It can be seen from these figures, as well as from the 

previous argument, that as z max becomes very large, these spectra take on a 

power-law form. Indeed, the spectrum predicted for the low-density model is 

compatable with the observed isotropic X-ray spectrum (Metzger, Anderson, 

van Dilla and Arnold 1964)"^ as well as a possible observation of the cosmic 

12 

gamma-ray flux above 100 MeV (Clark, Garmire and Kraushaar 1968), pro- 
vided the mean product of interacting matter and antimatter densities (n n-” 4 ) 
is of the order of 10” 18 cm -6 .* 


*Gamma-rays resulting from metagalactic cosmic -ray collisions will produce spectra similar to an- 
nihilation spectra, which peak near 70 |l + z m ax) 1 MeV as will be discussed elsewhere. The 
most noticeable observational difference between these spectra lies in the region above several 
hundred MeV where the annihilation spectrum drops sharply to zero while the collisional spectrum 
will have a power-law form. It may also be shown that the metagalactic gamma -ray spectrum re- 
sulting from cosmic -ray collisions will peak above or near 1 MeV and therefore cannot account for 
the observed isotropic X-ray spectrum. 
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The explanation most often considered for the isotropic X-radiation, is 

that it arises through Compton interactions of high-energy electrons and blackbody 

13 

photons in the metagalaxy (Felten and Morrison 1966). However, such an ex- 
planation is beset by various difficulties, as has been previously pointed out 

14 IK 

(Fazio, Stecker and Wright 1966; Brecher and Morrison 1967 ). 

The author would like to gratefully acknowledge Mr. Joseph Bredekamp for 
programming some of the numerical calculations presented in this paper. 
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FIGURE CAPTIONS 


Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


The normalized local differential gamma-ray spectrum from 
proton -antiproton annihilation. 

The metagalactic differential annihilation spectrum given for various 
maximum redshifts for an Einstein-de Sitter universe as discussed 
in the text. 

The metagalactic integral annihilation spectrum given for various 
maximum redshifts for an Einstein-de Sitter universe as discussed 
in the text. 

The metagalactic differential annihilation spectrum given for various 
maximum redshifts for a low density universe as discussed in the 
text. 

The metagalactic integral annihilation spectrum given for various 
maximum redshifts for a low density universe as discussed in the 


text. 



f r (E r ) 



Figure 1. The normalized local differential 
gamma-ray spectrum from proton-antiproton annihilation.. 
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